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Abstract 

The dissociative adsorption of methane on variously oxidized Pd, Pt and Pd-Pt 
surfaces is investigated using density-functional theory, as a step towards under- 
standing the combustion of methane on these materials. For Pd-Pt alloys, models of 
surface oxide structures are built on the basis of known oxides on Pd and Pt. The 
methane adsorption energy presents large variations depending on the oxide struc- 
ture and composition. Adsorption is endothermic on the bare Pd(lll) metal surface 
as well as on stable thin layer oxide structures such as the (y/E x y/E) surface oxide on 
Pd(100) and the Pt02-like oxide on Pt(lll). Instead, large adsorption energies are 
obtained for the (100) surface of bulk PdO, for metastable mixed Pdi-zPt^O^ oxide 
layers on Pt(100), and for Pd-Pt(lll) surfaces covered with one oxygen monolayer. 
In the latter case, we find a net thermodynamic preference for a direct conversion of 
methane to methanol, which remains adsorbed on the oxidized metal substrates via 
weak hydrogen-bond interactions. 

To whom correspondence should be addressed: arezoo.dianat@nano.tu-dresden.de 
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1 Introduction 



The catalytic oxidation of hydrocarbons in general, and of methane in particular, is con- 
sidered as an effective method of power generation associated with low emissions of CO2 
and N0 X . Because of the high H/C ratio of methane, the heat of combustion per mole 
of generated CO2 is higher than for other fuels, e.g. twice as much as for coal pQ. The 
interest towards the development and optimization of novel catalysts for the combustion of 
methane has thus considerably increased over the last years. Especially palladium-based 
catalysts have been extensively explored because of their high catalytic oxidation activity. 
Among those catalysts, the bimetallic system Pd-Pt has been the object of many investi- 
gations [2TTTT] . In several studies it was found that Pd-Pt catalysts exhibit higher methane 
conversion efficiency [5HB1I8HT2~| and better long-term stability [IIJ[7J[TT] than pure Pd. 

The performance of Pd-Pt catalysts strongly depends on the chemical state of the metal 
surface at the conditions suitable for oxidation, which in many cases lead to the formation 
of superficial oxide phases. For pure Pd and Pt catalysts, there is agreement that oxide 
formation has a positive effect on their methane oxidation activity [T5|[H]. although the 
active oxide phases have not been unequivocally determined yet. On the one hand, some 
authors have proposed that bulk PdO is less active than either a thin oxide layer [15] or a 
layer of adsorbed oxygen on Pd P^, [T7]. On the other hand, the formation of bulk oxide 
has been suggested to explain the observed increase in catalytic activity [THEE]. Recently, 
Gabasch et al. [20] reported that bulk PdO seeds grown on a surface, otherwise covered 
by a Pds04 surface oxide, are the active phase for methane oxidation on a Pd(lll) single 
crystal surface. The catalytic behavior of pure Pt metal is quite different. According 
to investigations in Refs. [T8l [19], [2TJ [22] , its catalytic activity reaches a maximum at a 
submonolayer coverage of adsorbed oxygen, and further oxidation with formation of Pt02 
leads to activity loss. Still, it is unclear whether a similar behavior can be expected also for 
other Pt oxide phases that could develop in an oxidizing atmosphere. Platinum bulk oxide 
phases comprise a-Pt02, /3-Pt02, Pts04 and PtO [23H25] . Correspondingly, a large variety 
of surface oxides can be expected [2"oT - |2"%] . Previous calculations have suggested that Pt 3 4 
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might be an active phase for the catalytic oxidation of carbon compounds [21]. Similarly, 
in the case of palladium, besides PdO also other bulk oxide phases could form, as e.g. 
Pd0 2 [29]. Different surface oxides have been observed on Pd(lll) [3QH32], Pd(100) [33] 
and on stepped surfaces [34] . 

Compared with the large amount of data available for pure Pd and Pt, little is known 
about the oxidation behavior of the bimetal Pd-Pt and its relation to the higher catalytic 
activity displayed by these alloys. Persson et al. [10] suggested that Pd-Pt catalysts on 
alumina consist of a PdO-rich phase coexisting with a palladium-rich Pd-Pt alloy. Studying 
the composition dependence of the catalytic activity of Pd-Pt, Lapisardi et al. [TT] found 
the highest activity for catalysts with very high Pd content Pdo.93Pto.07/Al2O3. 

In a previous study [35], we have investigated by means of density-functional theory 
(DFT) the thermodynamic stability of mixed bulk oxides Pdi-^Pt^Oy isostructural to 
known bulk oxide phases of pure Pd and Pt. According to these calculations, the mixed 
oxide phase Pdi_ x Pt x 02 with crystal structure analogous to a-Pt02 is stable only at 
rather low temperature (< 500 K for atmospheric oxygen pressure). With increasing 
temperature, a phase mixture of PdO with first Pt02 and later Pts04 becomes stable. At 
higher temperature, the oxides decompose directly to metallic Pd-Pt, except for Pd-rich 
systems, where PdO and metallic Pt coexist within a small temperature window. 

In the present DFT study, in a first step we analyze the stability of oxidized Pd-Pt 
surfaces as a basis for investigating their interaction with methane. Because of the lack 
of relevant experimental information, known surface oxide structures of pure Pd and pure 
Pt are used as guidelines to construct Pd-Pt surface oxide models. In a second step, we 
compute the driving forces for the dissociative adsorption of CH4 on various oxidized Pd-Pt 
surfaces, which is commonly believed to be the rate-determining reaction in the catalytic 
combustion of methane. 

The paper is outlined as follows. Computational details of the DFT analysis are briefly 
described in Sect. 2. In Sect. 3, we report on the calculated oxygen binding energies of 
various Pd-Pt oxide structures. The calculated methane adsorption energies on these oxide 
structures are then presented in Sect. 4. Finally, our results are discussed and interpreted 
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in relation to experimental findings in Sect. 5. 

2 Computational details 

Our DFT calculations are performed by means of the Vienna ab initio simulation pack- 
age (VASP) [5BH3E], using the PBE generalized gradient approximation (GGA) for the 
exchange-correlation functional [39J and the PAW method [HUE] • The wave functions are 
expanded in plane waves up to a kinetic energy cut-off of 400 eV. The periodically repeated 
simulation cells include slabs of six substrate layers covered with either adsorbed oxygen 
or a thin oxide layer, and with adsorbed CH 3 and H. In all simulations, the vacuum gap 
between the slab surface models is larger than 15 A. Unless stated otherwise, the size of 
the simulation cell corresponds to a (2x2) surface unit cell of the metal substrate. Inte- 
gration in the first Brillouin zone is performed using Monkhorst-Pack grids [12] including 
25 /c-points in the irreducible wedge. In all calculations, the positions of all atoms are 
optimized until all force components become less than 0.01 eV/A. Convergence of energy 
differences with respect to the used cut-off energies and Appoint grids is ensured within a 
tolerance of 10 meV/atom. Further computational details can be found in Ref. |43j . 

3 Surface oxide structures on Pd-Pt 

The oxidation of Pd and Pt surfaces proceeds from the chemisorption of oxygen atoms 
through the formation of surface oxides to the development of bulk oxide. In the case 
of palladium, the structure of surface oxide phases has been the subject of many experi- 
mental jUHlT] as well as theoretical [30 | l33 | l34 ^ l43 |l4^1 - l50] investigations. Phase diagrams 
of surface oxide structures in dependence on the chemical potential of oxygen have been 
thoroughly characterized for several surface orientations [3"0"|I3"3]. On Pd(lll), chemisorbed 
oxygen at low coverage is arranged with a p(2x2) periodicity. At higher coverage, several 
surface oxide phases form and coexist, as observed in STM investigations accompanied by 
theoretical modelling [S01S31SH] - On the contrary, on Pd(100) only one surface oxide has 
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been found [33JII5], consisting of two unit cells of PdO(lOl) over a (V5 x y/E) Pd(100) 
cell. Investigations of oxide formation on platinum revealed the following. On Pt(lll), 
the formation of a bulk-like, strongly distorted a-PtC>2 surface oxide was observed at an 
oxygen partial pressure of 0.5 atm and temperatures from 520-910 K [5HE2]. a-PtC^ was 
predicted to be the stable low temperature phase also by DFT calculations [23ES]- On 
Pt(100), DFT calculations suggest the formation of an a-Pt02-like surface oxide at low 
temperature [27l|28] and the existence of a stability region for a Pts04-like oxide layer at 
higher temperature. 



3.1 On-surface and sub-surface adsorbed oxygen 

To characterize the binding strength of oxygen atoms adsorbed at metal surfaces, we 
calculate an average oxygen binding energy Ef, per O atom according to the formula 



Eb = ~k 



Eq@s — E s 2^°2 



(1) 



where E Q @s is the total energy of the oxygen-metal system, Es the energy of the bare 
metal substrate, Eq 2 the energy of a free oxygen molecule, and No the number of O 
atoms. The oxygen atoms in formula ([I]) include oxygen located on the surface as well 
as in sub-surface positions. The calculated oxygen binding energies on (111) surfaces are 
listed in Tabled] for different oxygen coverages up to 1 monolayer (ML). The values in the 
last column are obtained for a stack of atomic layers with Pd in the surface layer, Pt in 
the sub-surface layer, and Pd in the remaining layers. In a previous DFT study [35], we 
have calculated oxygen adsorption energies on mixed Pd-Pt surface and sub-surface layers. 
For the considered Pd-Pt compositions and configurations, the oxygen adsorption energies 
have been found to vary between -1.26 eV and -0.94 eV. In the case of stacks of pure Pd 
and Pt layers, we have obtained an adsorption energy of -0.99 eV on Pt/Pd/Pd(lll) and 
-1.26 eV on Pd/Pt/Pd(lll). The latter value presents the strongest oxygen binding that 
we have found on all considered stacks of mixed metal layers. This strong binding has 
been attributed to charge transfer from Pd to Pt which causes a strong binding between 
the more electropositive Pd atoms and electronegative O atoms. 



The most stable oxygen adsorption sites on the (111) surface of Pd and Pt are the fee 
hollow sites [q*3"|I5"1]. According to the energy values in Table 1, the oxygen binding energies 
on pure Pt(lll) are significantly smaller than on pure Pd(lll) for all considered coverages. 
The strongest oxygen binding is found for the Pd/Pt/Pd(lll) layer stack. 

Further calculations are performed for the case of fixed sub-surface oxygen coverage of 
0.25 ML, and increasing on-surface coverage up to 0.75 ML, thus giving a maximum total 
coverage of 1 ML. Between the first and second metal layer there are three high-symmetry 
positions: an octahedral site underneath the fee on-surface hollow site, a tetrahedral site 
(tetra I) below the hep on-surface hollow site, and a second tetrahedral site (tetra II) 
directly below a surface metal atom [53]. For simultaneous on-surface and sub-surface 
adsorption, the fee hollow site on the surface and the tetra I sub-surface position are 
found to be the most stable oxygen positions for all considered oxygen coverages. The 
corresponding average oxygen binding energies are given in Tabled] (referred to as on+sub). 
They reveal that, in addition to on-surface adsorption, oxygen incorporation underneath 
the surface layer becomes favorable for a total oxygen coverage 9 tot > 0.75 ML on Pd(lll), 
and for 9 tot = 1 ML on Pt(lll). This result is in agreement with other DFT calculations 
for Pd(lll) [53]. For solely on-surface adsorption, the binding of oxygen is slightly stronger 
on Pd/Pt/Pd(lll) than on Pd(lll), whereas with sub-surface oxygen the average binding 
energy is larger for pure Pd(lll) due to the unfavorable binding of oxygen to the Pt 
sub-surface layer in the case of Pd/Pt/Pd(lll). 

3.2 Thin oxide layers 

With increasing oxygen coverage, surface oxide phases start forming on Pd and Pt surfaces. 
To our knowledge, there is no experimental information concerning the structure of such 
phases on the Pd-Pt bimetal. For this reason, model structures of possible Pd-Pt surface 
oxides are constructed here on the basis of the known oxide structures of pure Pd and 
Pt. In particular, we consider the experimentally observed PdO(101)-like [331115] and the 
theoretically predicted Pt 3 4 -like [27] oxide layers on the (100) surface, as well as the a- 
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PtCVlike [271111] layer on the (111) surface. According to our previous calculations 
these layers are characterized by small lattice misfits between oxide and metallic substrate, 
namely 2% for PdO(101)/Pd(100), 3.5% for a-PtO 2 (0001)/Pt(lll), and less than 1% for 
Pt3O4(100)/Pt(100). The corresponding oxide structures obtained after DFT geometry 
optimization are shown in Fig. [TJ 

In the case of PdO(101)/Pd(100), the Pd atoms of the oxide are located above the sur- 
face hollow sites of the substrate. The oxygen amount in the oxide layer corresponds to a 
coverage of 0.8 ML, referred to the surface metal atoms. In the case of a-PtO2(0001) /Pt(lll) 
the Pt atoms of the oxide are located above top sites, fee hollow sites, as well as hep hollow 
sites of the substrate, and the oxygen coverage is 1.5 ML. The epitaxial PtsO^lOO) /Pt(100) 
oxide layer corresponds to an oxygen coverage of 1 ML. In this case, the oxide layer includes 
two metal layers (cf. Figs. [T](e) and (f)). 

In the case of mixed Pd-Pt systems, we consider special Pd-Pt configurations in the 
oxide layers and in the metal surface layers which are generated by the following replace- 
ments. For example, starting from a pure Pd system in the case of Fig. la, we first replace 
Pd atom 1 in the oxide layer by Pt. With increasing Pt concentration in this layer Pd 
atoms 2, 3 and 4 are replaced by Pt. The Pt concentration in the metal surface layer is 
increased in the same manner starting with Pd atom 5. For the other two layer systems 
(Figs, lc and le), the Pt atoms are replaced by Pd in analogous manner. 

In order to characterize the stability of oxide structures with different Pd-Pt config- 
urations and variable oxygen coverage, an average oxygen binding energy can be defined 
as 

i r ... , ... at^ i 

(2) 



No 



T? TP ArLayer rpbulk AjLayer rpbulk No rp 
&Layer@S ~ ~ iV Pd ^Pd ~ iV Pi ^Pt 2~ 



where Eiayerms is the total energy of the substrate together with the oxide layer, Es the 
energy of the substrate alone, Epj and E b pl k the energies of Pd and Pt atoms in their bulk 
phases, and No, Np^ yer and Np^ yer the numbers of oxygen, Pd, and Pt atoms in the oxide 
layer. The binding energy (j2J) corresponds to the energy gained by deposition of molecular 
oxygen and formation of an oxide with metal atoms from ideal bulk metal reservoirs. With 
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this definition we neglect the difference of chemical potentials of the metal atoms in the 
pure metal bulk and in a Pd-Pt alloy. Since the free energy of mixing of Pd-Pt is of the 
order of few tens of meV per metal atom, while the oxygen binding energies are of the 
order of one eV, this approximation will not significantly influence our conclusions. 

To study the effect of the Pd-Pt composition on the oxygen binding energy, we vary 
the Pd/Pt ratio both in the oxide layer and in the outermost substrate- layer. All other 
substrate layers are composed of Pd in the case of the PdO(lOl) layer on Pd(100), and of Pt 
for the other two cases. The calculated oxygen binding energies for the three oxide layers are 
shown in Fig. [2] as a function of the Pt concentration of the outermost substrate layer. The 
different curves presented correspond to different compositions of the oxide layer. Every 
curve point represents one special Pd-Pt configuration as noted above. Test calculations for 
another configuration at fixed composition for Pdo.5Pto.5O on Pdo.4Pto.6/Pd(100) and for 
a-Pdo.33Pto.67O2 on Pdo.75Pto.25/Pt(lll) yield only small differences in the average oxygen 
binding energy of 0.012 eV and 0.006 eV, respectively. 

For the PdO(lOl) layer structure (Fig. 12(a)), the maximum variation of the oxygen 
binding energy with the substrate-layer composition is less than 0.05 eV. Similarly, the 
binding energy varies less than 0.1 eV for the a-Pt02 structure (Fig. [2(b), note the different 
energy scale). Replacement of Pd atoms of the PdO-like layer with Pt is energetically 
unfavorable. Namely, the O binding energy changes from -1.15 to -0.90 eV passing from 
a pure PdO to a pure PtO layer. The same holds for replacing Pt atoms of the a-Pt02 
layer with Pd. An analogous result has previously been found for mixed bulk oxides of 
Pd-Pt [35]. For the case of the Pt304-like oxide layer (Fig. 12(c)), the changes in the 
binding energy with the substrate-layer composition are also small (less than 0.12 eV). 
For this oxide layer it is very unfavorable to replace Pt atoms with Pd. Indeed, a Bader 
analysis [551156] reveals that the oxygen atoms in Pd3O4/Pt(100) gain 0.67 e from metal 
atoms, while the corresponding value for PtsO4/Pt(100) is 0.73 e, reflecting the stronger 
oxygen binding in the latter case. For the other two oxide structures considered, the charge 
transfer from metal atoms to oxygen does not change significantly with variation of the 
Pd-Pt composition in the oxide layer, consistently with the calculated small differences in 
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the oxygen binding energy. 



4 Dissociative adsorption of methane 

The first important step in the catalytic oxidation of methane is its adsorption on the 
catalyst surface and dissociation into adsorbed methyl and hydrogen. To get a first in- 
sight into the catalytic activity of Pd-Pt catalysts for methane oxidation, we compute 
the thermodynamic driving force for this adsorption reaction on various oxide structures, 
in particular on the superficial oxides of Pd-Pt considered in the previous section. The 
adsorption energy for the dissociative adsorption of CH 4 to CH 3 and H is defined as 

E a = EcH 3 k,H®s — Es — Ech a , (3) 

where Ech 3 &h@s is the total energy of methyl and hydrogen adsorbed on the substrate, 
Es the energy of the substrate, and Ech a the energy of methane in the gas phase. We do 
not perform zero point energy corrections of the calculated adsorption energies since we 
think that these corrections lead to comparatively small shifts of energy values (see e.g. 
Ref. [57]), which is of minor importance for comparing adsorption energies on different 
systems as main goal of this study. 

4.1 CH4 adsorption on clean metal surfaces 

We first report on our calculations concerning methane adsorption on pure metallic (111) 
and (100) surfaces. In agreement with other DFT calculations [581159"] . we find that on the 
(111) surface the most stable adsorption site for CH 3 is on top of metal atoms and for H on 
fee hollow sites. This applies to pure Pd(lll) and Pt(lll) as well as to the Pd/Pt/Pd(lll) 
layer stack with Pd in the surface and Pt in the sub-surface atomic layer. The calculated 
adsorption energies, shown in Tabled are -0.08 eV for the Pt(lll) surface and 0.19 eV 
for the Pd(lll) surface. The positive value for Pd(lll) corresponds to an endothermic 
reaction. This is in agreement with previous DFT calculations for methane adsorption on 
Pd(lll) within the generalized gradient approximation (GGA) (E a = 0.27 eV), whereas 
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an exothermic reaction was found by using the local density approximation (LDA) (E a = 
-0.62 eV) [SH]- For comparison, at the LDA level we obtain adsorption energies of -0.43 eV, 
-0.66 eV, and -0.45 eV on Pd(lll), Pt(lll), and Pd/Pt/Pd(lll), respectively. Thus, for 
all considered systems, the LDA values are shifted roughly by -0.6 eV with respect to the 
GGA values. 

In the case of all (100) surfaces, CH3 on top of metal atoms and H on bridge sites are 
the most stable adsorption configurations. On Pd(100) and Pd/Pt/Pd(100), dissociative 
adsorption of methane is found to be endothermic, in agreement with DFT calculations 
in [5D]. Only on Pt(100), methane adsorption is exothermic with an adsorption energy of 
-0.31 eV. 

4.2 CH4 adsorption on surfaces with adsorbed oxygen 

After analyzing methane adsorption on clean metal surfaces, we investigate here the adsorp- 
tion on metal surfaces covered with 1 ML oxygen. A first set of calculations is performed 
starting with 1 ML of oxygen adsorbed solely on the (111) surface. As initial condition for 
the DFT calculations, the following positions of adsorbed CH3 and H are chosen: H always 
on top of oxygen atom, and (i) CH 3 on top of metal atom, (ii) CH 3 on top of oxygen atom, 
and (iii) CH 3 on hep hollow site. For these three cases, the calculated adsorption ener- 
gies after structural relaxation are listed in Table [2] The corresponding optimized atomic 
structures are shown in Figs. [3] and H] for the Pd(lll) and Pt(lll) surfaces. The structures 
obtained for the Pd/Pt/Pd(lll) layer stack are qualitatively the same as for Pd(lll). The 
largest adsorption energies are obtained starting with CH3 on top of a metal atom and 
H on top of an oxygen atom. In particular, the adsorption energy on Pd/Pt/Pd(lll) is 
slightly higher than on Pd(lll). In the latter two cases, structural optimization leads to a 
remarkable reconstruction of the adsorbate layer (Fig. 02(a)). Both the CH 3 group and the 
H atom move away from their initial adsorption sites and bind to the same O atom, thus 
forming an adsorbed methanol molecule. Notably, if CH3 and H are initially placed over 
other adsorption sites on Pd(lll), CH 3 and H also detach from the metal atoms, but bind 
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separately to different O atoms of the surface, resulting in adsorbed OH and CH 3 groups 
and no methanol formation. The same is found in the case of the Pt(lll) surface where 
CH3 and H remain bound to separate sites. In this case, the Pt atoms to which CH3 is 
bound are lifted off the surface layer by 2.2 to 2.5 A, depending on the initial configuration 
(cf. Fig.©. 

Analogous calculations of methane adsorption energies have been performed for surfaces 
covered with 0.75 ML oxygen on the surface and 0.25 ML in sub-surface positions. Again, 
three different initial positions of adsorbed CH 3 and H have been chosen: H always on top 
of oxygen, and (i) CH 3 on top of metal atom, (ii) CH 3 on top of oxygen atom, and (hi) 
CH 3 on fee hollow site (oxygen vacancy). The largest adsorption energy results for CH 3 
and H initially on top of oxygen for all considered surfaces (denoted by CH 3 -fcc & H-fcc in 
Table [2]). The energy values in Table |5] indicate that methane adsorption on surfaces with 
sub-surface oxygen is about 1 eV weaker than for the case of oxygen adsorbed solely on 
the surface. In the presence of sub-surface oxygen, methanol does not form spontaneously 
starting from the initial geometries above, H and CH 3 remaining separately bound to 
different oxygen atoms. 

While in the simulations so far the formation of methanol takes place only on Pd(lll) 
and Pd/Pt/Pd(lll) with oxygen on-surface coverage, we now calculate the adsorption 
energy after methanol formation on all other surfaces, irrespective of the initial geometry 
choosen. To this end, a novel set of DFT structural relaxations is performed starting 
with methanol adsorbed on all surfaces, including those with subsurface oxygen, using the 
atomic configuration obtained for Pd(lll). In all cases, the adsorption energies are higher 
than for separately adsorbed CH 3 and H, as reported in Table [2J As visible in Fig. |3fa), 
the methanol molecule remains loosely adsorbed to the surface, mainly via a hydrogen- 
bond between the OH group of methanol and a surface O atom. The methanol adsorption 
energies amount to -0.11 eV on Pd(lll) and Pd/Pt/Pd(lll), and -0.08 eV on Pt(lll) for 
1 ML on-surface oxygen coverage, reflecting weak molecule-surface interactions. 
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4.3 CH4 adsorption on superficial oxide layers 

With increasing oxygen supply, the oxidation of the metallic catalyst is expected to pro- 
ceed with formation of thin oxide layers. Correspondingly, we further investigate methane 
adsorption on the thin oxide layer structures described in Sect. 3.2. The most favorable 
adsorption sites of methyl and hydrogen are determined via relaxation of different struc- 
tures with the following initial positions of the adsorbates: both CH3 and H on top of 
oxygen or metal atoms, CH 3 on top of oxygen and H on top of metal atom, and vice versa. 
The optimized atomic structures on the PdO(101)/Pd(100), a-PtO 2 (0001)/Pt(lll), and 
Pt3O4(100)/Pt(100) layers are shown in Fig. [5j The Pd-Pt composition is varied both 
in the oxide layer and the outermost layer of the metal substrate (cf. Fig. [6]). For all 
compositions, the stable atomic configurations of adsorbed CH 3 and H have been found to 
be qualitatively equal. However, the values of the adsorption energy, presented in Fig. El 
show comparatively large variations. Every curve point in Fig. 6 represents a special Pd-Pt 
configuration in the oxide and metal surface layer as described above. Test calculations of 
another Pd-Pt configuration for a PdO(lOl)- and PtCvlike oxide layer at fixed composi- 
tion (see Sect. 3.2) yield only small differences in the methane adsorption energy of less 
than 0.05 eV. 

For the PdO(101)-like oxide layers, the most stable adsorption sites for CH 3 are ei- 
ther on Pd, if only Pd atoms are present, or on Pt, if Pt replaces Pd atoms in the oxide 
(Fig. E(a)). The corresponding adsorption energy values (Fig. EJ^a)) indicate that the 
methane decomposition reaction on this oxide structure is endothermic for all Pd-Pt com- 
positions considered. As a general trend, for a given composition of the oxide layer, the 
reaction is the more endothermic the more Pd atoms are present in the outermost layer of 
the metal substrate, except for pure PtO(lOl) where the adsorption energy is independent 
of the substrate composition. For a given substrate composition, increasing the Pd/Pt 
ratio in the oxide from to 1 results in a rather complex behavior of the adsorption energy 
values, whereby pure PtO(lOl) layer is the least reactive structure in all cases. 

In the case of the a-PtCvlike layer, where all metal atoms are fully coordinated by 
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oxygen, adsorption of both CH 3 and H occurs necessarily on the O atoms (Fig. 0(b)). 
The adsorption energy depends only weakly on the Pt concentration in the outermost 
substrate-layer, whereas a strong dependence on the oxide layer composition is observed 
(Fig. E(b)). While the reaction is clearly endothermic for the pure PtC>2 oxide layer, it 
becomes exothermic for a Pd/Pt ratio in the oxide larger than about 0.5. Thus, the driving 
force for methane dissociation increases monotonously with the Pd content in the oxide. 

In the case of the Pt304-like oxide layers, both CH3 and H adsorb strongly on the un- 
dercoordinated exposed oxygen atoms (Fig.H^c)), as found also previously for a pure Pt 3 04 
layer [27J. The adsorption energy of -1.74 eV obtained here is comparable to the value of 
-1.47 eV computed in Ref. [27J using norm-conserving pseudopotentials. Increasing the 
Pd concentration in the oxide layer leads to stronger and stronger adsorption, with little 
dependence on the composition of the outermost substrate-layer (Fig. E(c)). In particular, 
the large adsorption energy of -2.5 eV, computed for the PdsO4/Pt(100) oxide, is compa- 
rable to the driving force for methane dissociation and methanol formation computed for 
the case of an oxygen ML adsorbed on the Pd/Pt/Pd(lll) surface (see Table 2). 

4.4 CH4 adsorption on bulk PdO-like oxides 

After considering methane adsorption on thin oxide layers on Pd-Pt metal substrates, it 
is interesting to compute, for comparison, the methane adsorption energies on the corre- 
sponding bulk oxide surfaces. We note that for a-Pt02 the adsorption on a single layer 
is already representative of the behavior of the bulk oxide, which consists of stacked PtC>2 
layers weakly bound by van der Waals forces [27]. Furthermore, the surface of bulk Pta04 
with lowest surface energy is nearly identical to the surface of the thin oxide layer consid- 
ered above. In the following, we thus consider only surfaces of mixed Pd-Pt oxides with 
the structure of the well-known PdO bulk oxide phase. In particular, methane adsorp- 
tion energies are calculated on the (100) and (101) surfaces, which present low surface 
energies [50J. 

As a model of a mixed oxide Pd-Pt-O, we analyze a slab of five PdO layers covered with a 
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mixed Pd 1 „ ;c Pt x O layer, in which Pd and Pt are arranged in a chess-like pattern for x = 0.5 
(adsorption energies for the row-like Pd-Pt pattern differ less than 0.025 eV from the chess- 
like one). Our analysis shows that on the (100) surface CH3 and H adsorb preferentially on 
top of oxygen atoms, in agreement with previous calculations [61] (Fig. 0(a)). On the (101) 
surface, H adsorbs on top of oxygen and CH 3 on top of metal atoms, preferentially on Pt in 
the case of a mixed Pd-Pt bulk oxide, as also found for the thin oxide layer in the previous 
section (Fig. [3(b)). The calculated adsorption energies on the different surfaces are listed 
in Table [31 With increasing Pt content in the oxide, a monotonously decreasing driving 
force for methane adsorption is found on PdO(100), whereas on Pd(101) the adsorption 
energy displays a non-trivial dependence on the Pt content, similarly as in the case of the 
thin Pd(101)-like oxide layer considered above. Unlike the case of thin PdO-like oxide 
layers, where methane adsorption is endothermic, for the PdO bulk phase, the adsorption 
reaction is exothermic. Regarding the dependence of the adsorption energy on the oxide 
composition, the largest value of -1.0 eV is obtained for the (100) surface of pure PdO (cf. 
Fig. G3(a)). This value is however considerably lower than the energy values calculated for 
the Pt 3 4 -like oxide layers on Pt(100) and for the oxygen adlayers on Pd(lll) (cf. Fig.^c) 
and Table 2). 

5 Discussion 

The dissociative adsorption of methane on noble metal surfaces represents an essential 
reaction step of the catalytic combustion of methane, and is thought to limit the reaction 
kinetics. As a first effort towards a deeper understanding of this catalytic reaction, we have 
studied methane adsorption on various oxide structures of Pd, Pt and Pd-Pt alloy surfaces. 
The choice of the investigated systems is motivated by the facts that palladium oxide 
phases, in particular PdO, have been suggested to be catalytically more active than pure 
Pd [IHMIH], and that mixed Pd-Pt catalysts have been found to possess higher conversion 
efficiency for methane combustion, especially on the long term. Our investigation includes 
chemisorbed O atoms, thin oxide layers and bulk PdO-like oxide surfaces, addressing both 
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their thermodynamic stability and their reactivity towards methane dissociation. 

5.1 Formation of oxide layers on Pd-Pt surfaces 

The thermodynamic stability of different oxidized Pd-Pt structures has been addressed 
by calculating average oxygen binding energies at K. In general, superficial oxide layers 
are found to be more stable than oxygen adlayers(including sub-surface O atoms) at the 
corresponding oxygen coverages, as displayed in Table 2 and Fig. [2J The larger stability 
of oxygen adlayers on Pd compared to those on Pt can be understood simply from the 
lower electronegativity of palladium. Consistently with previous findings [13], the largest 
binding energy is found for a Pd/Pt/Pd(lll) layer stack, due to the partial donation of 
electrons from the Pd surface layer to the Pt sub-surface layer, which increases the surface 
reactivity. 

In the case of surface oxide layers, the Pd-Pt composition of the outermost layer of 
the metal substrate has only a minor influence on the oxygen binding strength. Changes 
of the binding strength due to varying the Pd/Pt ratio in the oxide can be explained on 
the basis of the stability of the corresponding bulk oxide phases, as thoroughly addressed 
in Refs. [27JE5]- Namely, mixed oxides with PdO structure are destabilized by replacing 
Pd atoms with Pt, while mixed oxides with a-Pt02 or Pts04 structure are destabilized by 
replacing Pt atoms with Pd. 

5.2 Dissociative adsorption of CH4 on oxidized Pd-Pt surfaces 

Concerning the dissociative adsorption of methane, our DFT calculations at the GGA 
level suggest, in agreement with existing literature [5H], that Pd metal surfaces are not 
reactive, and Pt surfaces are only little reactive. Since CH3 and H possess an electron- 
donor character, they are expected to bind better to electronegative elements. Indeed, 
binding to pure Pt (electronegativity 2.3 on the Pauling scale) is stronger than to pure Pd 
(electronegativity 2.2). However, if adsorbed oxygen (electronegativity 3.4) is present on 
the surface, either in the form of an oxygen adlayer or as surface oxide, then binding to 
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oxygen provides a strong driving force for methane dissociation. 

Correspondingly, we have obtained the largest adsorption energies for 1 ML of oxygen 
adsorbed on the (111) metal surfaces, in particular on Pd(lll) and on the Pd/Pt/Pd(lll) 
layer stack, with an energy gain of about -2.4 eV (cf. Table 2). A peculiar effect observed on 
the latter two surfaces is the spontaneous formation of a methanol molecule which remains 
only loosely bound to the surface via hydrogen bonds and weak metal-oxygen interactions 
(Fig. |3£a)). We observe a net thermodynamic preference for methanol formation on all 
three surfaces considered, with the largest values being obtained on Pd(lll). An analysis 
of the energy barriers associated with the CH 4 — > CH 3 OH conversion reaction exceeds the 
scope of the present investigation. However, the possibility of a direct methane to methanol 
conversion on the (111) surface of transition metals has been recently put forward in a 
theoretical DFT study [63J. It has to be noted, however, that the further combustion 
of methanol on Pd-Pt surfaces is expected to take place at the temperatures required 
to dissociate the C-H bonds of methane. Therefore, in general the selectivity towards 
methanol formation is found to decrease by increasing the overall conversion efficiency 

In the case of superficial oxide layers formed on Pd-Pt, we have found that methane 
adsorption is clearly endothermic both on the thin PdO(101)-like layer on (100) surfaces 
and on the thin a-PtCVlike layer on (111) surfaces. Since the same was noted for oxidation 
of carbon monoxide [22], these phases may be considered as a kind of passivation layer on 
the metal substrate, whose formation may suddenly reduce the oxidation activity of the 
catalysts. This property, however, is not shared by the surfaces of bulk oxides, where both 
the (101) and the (100) surfaces are reactive towards methane dissociation (Table 3). This 
behavior can be understood by considering the charge transfer between the CH 3 molecule 
and the surface atoms. According to a Bader analysis, in the case of bulk PdO, the CH3 
molecule donates electrons to the surface and becomes positively charged (+0.44e), while 
in the case of the PdO-like oxide layer CH 3 gains electrons (-0.13e). This indicates a 
strongly reduced electronegativity of the metal atoms of the thin oxide layer due to the 
presence of the underlying metal substrate, compared with the surfaces of bulk oxides. As 
a consequence, binding of CH 3 to the surface of the bulk oxide is favorable, in contrast to 
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the oxide monolayer over the metal substrate. This finding is consistent with experimental 
observations that oxidation of Pd with formation of PdO is beneficial for methane oxidation 
[18l[T9]. 111 comparison, bulk a-Pt02 is inert with respect to methane dissociation in the 
absence of defects (see also Ref. [2Z])- This again is consistent with observations that 
formation of bulk oxide phases is beneficial for methane oxidation on Pd, but not on 
Pt [TB1HH1EI]- Oxidation of Pt might be beneficial if the reaction conditions allow the 
formation of Pt 3 04-like phases (see also the discussion in [2U|27J|28]). Namely, our DFT 
calculations predict adsorption energies of CH 3 and H on Pt 3 04 and Pd 3 04 which are 
higher than those on PdO surfaces and comparable with that on oxygen adlayers covering 
Pd(lll). This is due to the availability of undercoordinated oxygen sites on the exposed 
Pt 3 O4(100) surface, to which electron acception from CH 3 and H is very favorable. 

Concerning effects on methane adsorption resulting from Pd-Pt alloying, we note that in 
most of our calculations the composition of the metal substrate underneath the thin oxide 
layer does not show a pronounced effect on the methane adsorption energy. However, 
the effect of the composition of the oxide layer is evident and the reactivity decreases with 
increasing Pt content. An exception is the thin PdO(101)-like oxide layer on (100) surfaces, 
where the CH 3 group binds preferentially to the metal atom rather than to oxygen. In 
this case, the presence of Pt atoms does increase the reactivity since they are stronger 
electron acceptors, as mentioned above. However, at the same time the presence of Pt 
results in a decrease of the reactivity of the O atoms towards the adsorption of H. These 
two counteracting effects result in a non-trivial trend of the computed adsorption energies 
with increasing Pt content, both in the case of the thin layers (Fig.E]) and of bulk PdO(lOl) 
(Table 3). The actual adsorption energies in this case depend on the specific arrangements 
of Pt, Pd and O atoms close to the adsorption sites of CH 3 and H, both in the oxide layer 
and in the metal substrate. However, in none of the cases considered is the adsorption 
reaction exothermic, as discussed previously. 
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6 Conclusions 



In conclusion, we have performed extensive calculations of the driving force for methane 
dissociation on clean and oxidized Pd-Pt surfaces. Both in the case of Pd and Pt, formation 
of stable thin layer oxide structures such as the PdO(101)-like layer on the (100) surface and 
the «-Pt02-like layer on the (111) surface, leads to a reactivity loss towards the dissociative 
adsorption of methane. Methane adsorption is instead favored on metastable surface oxide 
structures such as a PtsC^-like layer, whose reactivity increases with increasing Pd content. 
Furthermore, reactivity is recovered for exposed surfaces of bulk PdO, consistently with 
existing experimental results [6|[TTj. 

In the case of oxygen adlayers we have found that the formation of a methanol molecule 
after methane dissociation is thermodynamically favored. This may suggest that selective 
conversion of methane to methanol rather than total oxidation can be achieved under 
pressure and temperature conditions that prevent the formation of surface oxides. However, 
the further combustion of the produced methanol over other sites of the Pd-Pt surface 
cannot be excluded, reducing the selectivity of the conversion reaction. 

As a final remark, we would like to note that our investigation, being concerned with the 
thermodynamic driving force for the methane dissociation reaction, should be considered 
only as a preliminary indication about the reactivity of different surface structures. Based 
on the results presented here, more thorough calculations of the corresponding activation 
barriers, which can be directly related to the kinetic constants for adsorption, shall be the 
subject of forthcoming works. 
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Tables 



Table 1 : Calculated average oxygen binding energies (in eV) in the case of only on-surface 
adsorption (on) and for simultaneous on-surface and sub-surface adsorption (on+sub) on 
Pd(lll), Pt(lll), and Pd/Pt/Pd(lll) at different O coverages. In the on+sub case, the 
amount of sub-surface oxygen is always 0.25 ML and the remaining oxygen is on the surface, 
i.e. for 0.25 ML total coverage there is no oxygen on the surface in the on+sub case. 



Total coverage 
in ML 


Pd(lll) 
on on+sub 


Pt 
on 


111) 
on+sub 


Pd/Pt/Pd(lll) 
on on+sub 


0.25 


-1.34 


0.36 


-1.22 


0.68 


-1.39 


0.64 


0.50 


-1.03 


-0.93 


-0.92 


-0.52 


-1.06 


-0.85 


0.75 


-0.61 


-0.72 


-0.55 


-0.47 


-0.69 


-0.68 


1.00 


-0.22 


-0.42 


-0.17 


-0.27 


-0.25 


-0.41 
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Table 2: Calculated adsorption energies (in eV) for the dissociative adsorption of methane 
on clean (111) and (100) surfaces of Pd, Pt, and Pd/Pt/Pd, and on surfaces covered with 
1 ML oxygen on the surface, as well as with 0.75 ML on the surface and 0.25 ML in sub- 
surface positions. For clean surfaces, the values in squared brackets (in italics) correspond 
to adsorption energies within LDA. In the calculations, the initial positions of CH 3 and 
H have been chosen on high-symmetry sites of the (111) metal surface: CH3 on top of 
metal atom (-top) and H on fee hollow site on top of oxygen (-fee) in case of the clean 
surfaces. For the systems with sub-surface oxygen, two fee sites have been considered, CH 3 
on fee hollow site on oxygen atom (-fee) and CH3 on oxygen vacancy (-fee*). The values in 
parentheses are the adsorption energies of methanol formation on the surfaces, irrespective 
of the reaction path. (These value are not related to the initial CH 3 and H positions in 
the table, and are arbitrarily reported besides the corresponding highest energy values 
obtained for the other initial configurations). 



Clean surfaces 


Pd(lll) 


Pt(lll) 


Pd/Pt/Pd(lll) Pd(100) 


Pt(100) 


Pd/Pt/Pd(100) 


0.19 [-O.43J 


-0.08 [-0.66] 


0.22 [-O.45] 0.29 


-0.31 


0.35 


1 ML oxygen on (111) surfaces 



CH 3 -top & H-fcc CH 3 -fcc & H-fcc CH 3 -hcp & H-fcc 



-2.42 -1.65 -1.39 

2.13 (-2.27) -1.62 -1.99 

-2.47 -1.74 -1.70 

0.75 ML oxygen on-surface/0.25 ML oxygen sub-surface 

CH 3 -top & H-fcc CH 3 -fcc & H-fcc CH 3 -fcc* & H-fcc 

0.80 -1.59 (-1.87) -0.75 

1.39 -1.40 (-1.78) -1.23 

0.95 -1.28 (-1.88) -0.61 



Pd(lll) 
Pt(lll) 
Pd/Pt/Pd(lll) 



Pd(lll) 
Pt(lll) 
Pd/Pt/Pd(lll) 
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Table 3: Calculated methane adsorption energies (in eV) on the (101) and (100) surfaces 
of the PdO bulk oxide phase with one mixed oxide surface layer Pdi_ x Pta.O. 



Structure 


x=0 


x=0.25 


x=0.5 


x=0.75 


x=l 


PdO(lOl) 


-0.37 


-0.63 


-0.82 


-0.61 


-0.45 


PdO(lOO) 


-1.00 


-0.94 


-0.57 


-0.51 


-0.18 
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Figures 




Figure 1: Top- and side-views of optimized atomic structures of different oxide layers: (a,b) 
PdO(lOl) on Pd(100), (c,d) a-PtO 2 (0001) on Pt(lll), and (e,f) Pt 3 O 4 (100) on Pt(100) (O 
- red small spheres, metal - large spheres). In the case of mixed oxide and surface metal 
layers, Pd atoms are replaced by Pt (and vice versa) as described in the text. 
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Figure 2: Calculated average binding energies of oxygen for different oxide layer structures 
as a function of the Pt concentration x of the outermost substrate layer. The different 
curves correspond to different compositions of the oxide layer, as indicated by the labels, 
(a) Pdi_^Pt y O(101) onPdi_ x Pt a; /Pd(100), (b) a-Pdi_^Pt y O 2 (0001) on Pdi_ s Pt :c /Pt(lll), 
and (c) Pdi_ 2/ Ptj / O 4/3 (100) on Pd 1 _ x .Pt a ./Pt(100). 
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Figure 3: Optimized atomic structures of adsorbed CH 3 and H on a Pd(lll) metal surface 
covered with 1 ML oxygen. In the calculations the initial lateral positions of CH 3 and H 
have been chosen as follows: a) CH3 on top of a metal atom and H on top of oxygen, b) 
CH 3 and H on top of oxygen atoms, and c) CH 3 on hep hollow site and H on top of oxygen. 
The largest adsorption energy is obtained for a). Distances between adsorbed atoms and 
nearest neighbor substrate atoms in A: a) dc-pd = 4.00, dc-o = 1-43, dn-o = 0.98; b) 
d c ~Pd = 3.04, d C -o = 1-41, d H _o = 1-00; c) d C - Pd = 3.05, d c ~o = 1-40, d H _ Q = 0.98. 
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Figure 4: Optimized atomic structures of adsorbed CH 3 and H on a Pt(lll) metal surface 
covered with 1 ML oxygen. In the calculations the initial lateral positions of CH 3 and H 
have been chosen as follows: a) CH3 on top of a metal atom and H on top of oxygen, b) 
CH 3 and H on top of oxygen atoms, and c) CH 3 on hep hollow site and H on top of oxygen. 
The largest adsorption energy is obtained for a). Distances between adsorbed atoms and 
nearest neighbor substrate atoms in A: a) dc-pa = 2.97, dc-o = 1-42, dn-o = 0.99; b) 
dc-P d = 3.03, d C -o = 1-41, d H _o = 1-00; c) d C - Pd = 2.92, d C -o = 1-41, d H _ Q = 1-00. 
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Figure 5: Optimized atomic structures of adsorbed CH 3 and H on thin oxide layers of 
pure a) PdO(101)/Pd(100), b) a-PtO 2 (0001)/Pt(lll), and c) Pt 3 O 4 (100)/Pt(100). The 
shown atomic configurations correspond to the highest adsorption energy for each oxide 
layer. Distances between adsorbed atoms and nearest neighbor substrate atoms in A: a) 
d C - Pd = 2.05, d C -o = 2.77, d H _ Q = 0.98; b) d C - Pd = 3.15, d C -o = 1-43, d H . Q = 0.98; 
c) dc-P d = 2.85, d H _o = 0.98, d C -o = 1-44. 



31 




0.2 0.4 0.6 0.8 

Surface Pt concentration (x) 



PtCL 



b) 



bo 



Pd 0.33 Pt 0.66 O 2 



e 
_o 

& 
o 



Pd 0.66 Pt 0.33 O 2 



PdO, 



> -1.5 



0.25 0.5 0.75 

Surface Pt concentration (x) 



PtO_ 



c) 



Pd 0.33 Pt 0.67 O 4/3 



C 

O 



Pd 0.67 Pt 0.33 O 4/3 



-a 
< 



PdO, 



-2.5 - 



0.25 0.5 0.75 

Surface Pt concentration (x) 



Figure 6: Calculated methane adsorption energies on thin Pd-Pt oxide layers as a 
function of the Pt concentration x in the metal surface layer with the Pt concentra- 
tion y in the oxide as a parameter: (a) Pdi_j / Pt y O(101)/Pdi_ x Pt x ./Pd(100), (b) a- 
Pd 1 _ J ,Pt y 2 /Pd 1 _ x Pt x /Pt(lll), and (c) Pdi_ tf Pt w O 4 / 3 /Pdi_ x Pt x /Pt(100). 
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Figure 7: Optimized atomic structures of adsorbed CH3 and H on the a) (100) and b) 
(101) surface of the bulk oxide phase PdO. Distances between adsorbed atoms and nearest 
neighbor substrate atoms in A: a) dc-pa = 2.96, dc-o = 1-47, dn-o — 0.99; b) dc-pd = 
2.04, d c _o = 2.81, d H -o = 0.98. 
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